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Abst’^a jt 

The hunan operator can be de cribed by a linear nodel and a remnant 
added to the output of this model. In practical cases this remnant 
is always non-tero. Therefore the presence of this remnant has to 
te ta'Ken into account in any identification method applied in human 
operator re sear -h . 



Parameter estimation techniques are discussed with emphasis on un- 
bijscid estimates in the presence of nois-». A distinction between 
open and closed loop systems ic made. A method is given based on 
the application of extan.al forcing functions consisting of a sum 
uf sinusoids; this method is thus based on the estimation of Fourier 
coefficients and is applicable for mcdels with poles and zeros in 
open and closed loop systems. 


To describe the behavior of the human operator controlling time 
invariant linear system with time conetante up to the order of 
seconds, the Describing Function Method is a powerful too 
This DFM says that the system to be considered, in thie case the 
human operator's behavior, can be aseumed to eonelet of a linear 
time invariant system given by a transfer function and a remnant 
uncorrelated with the system input and added to the ouput of the 
system* The methods of estimating the parameters of these quasi- 
lin^^T mudels can he devlded it two main groups [2]* 

^ From the observation of input-ouput of a 

system, the Impulse response or t'.e transfer function is deter- 
mined. These methods deal with classical identification tech- 
niques, such as the determination of Bode or Nyquist plots from 
sinusoidal test signals. More reeeet additions are correlation 
techniques, power density analyses, FFT's and averaged response 
techniques* 

^ observing the input-output of a system, 

the parameters of a defined mathematical model can be estimated. 
Therefore both the system to be analysed and the mathematical 
model are given the same input. From the differences between the 
outputs of system and model the parameters can be determined 
according to a certain performance criterion and by means of a 
wall considered strategy. 

The funoamental difference between both the methods is that in a 
General Method only knowledge of the bandwidth of the input of the 
system to be a. ilyaed is required (to choose a proper test signal), 
while for a More Specific Hs'^hod the structure of the model 
itself must be Known as well. This pap« ^ will deal with the dis- 
cussion of a More Specific Method applied to the description of 
the human operator in open and closed loops. 
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As pointed out before, a iyaten in an epea loop ean bt daaaribed 
by e linear transfer function H(p) and a remnant n(t) unaarralatad 
ath the system input. For a known tranafar funetien H(p) tha 
parameters of a model baaed on this tranafar funetien ean be eaX* 
culated from a set of equations which reauita from minimiainp a 
function Ctc(t)it> [l], where tha quantity c(t) la the diffaranea 
between system output y(t) end model output y^^t) (aae Flfi !)• 



FiaURB u 

Blook diagram of ayatam idontifiaation bp maana Pf 
parameter pptimationi p demolaa thp Mpppipidp ppppptpr 

d/dt. 

If the function E(c(t>;t) ia equal to the moan aquarad error 

1 T 2 

T /(c(t)] dt then also the error c(t) will be unaarralatad with 

0 

the s> >t*em £nput x(t). In the apeeial oaaa that tha model output 
ean be written aes 

»•<*» * ih <»> 


SOS* 


• the oereeetepe e. «.n he eeieyleted free • eet el lieeee e«Mtle«et 
■edele heviht thle ppepepty ere eelXed *Xl»eep 1» the shpeootepe". 

A eyetee heving e trenefer fyoetlee with lepee eeljr ie •• e«.*epSe 
el e*eh a eedeH epiteet deeerihed h| e tranefep fnhetlM with 
pelee and neree era net linear in the parMetern. Hewewar. ler 
thaee ayatane a aa-*aallad "panarallnad nadal* (sj ean ha dafltad 
t— ri«. *). Hera, net tha dlffaranaa .<t) hatwaaa ayntan awtpnt yit» 
and Mddl antpwt y»it> la wialninad aaeardln* ta a siren orltarlan 
s(*(t)tt)« 



pMnatft. 


rmut »i 

Param$ur t$tima$itn Sp Mint a 
pandMllnaO naOdt> 

bat thn dlflarnnoa lit) In alnUlMd annardlng tn tho dbltablao 
SU(t)it)t th# swnntlty lit) In dnllnnd nt tb# dlffaonhnd botwam 
two alshdldi wblah are obtained by Mdllylnt tho tyataa iBpnt n(t> 
with tha trantfar Innatlan N(p) and tha awtpnt y(t) with tho tbbM* 

Par Iwnotlan Sip)i tha trahnfar Iwhotleht Rip) and Sip) pwbbotn 
only tarati Rawarnri It aan ht thawn that thd antlMtdb at tbo 
paraMtart ebtalntd in thld way will b# blaaod it tha error alsnal (it) 
ia nan«whlta tlj# tharalara, thla nathad aah ba appllod only II the 
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r«nn*nt nit) 1« t«ro or v«ry •b.aU in rolotlon to tht lyittfli «utyyt 
v(t). .f «n untinsoJ ontlnato It in moot prtttietl ««••• 

•<lth«r thtf ptratiottrt of a ttt of filttrt for tht tigntlt x(t) ang 
y( : i tava *^a cttiiaattd which whittn tht tiintl (it)» and thia hat 
ba tccoapithtd ittrativtly [3] » or tha nedaX paranatara hava to 
ba aatlnttad from tht original eritarion C(c<t>|t/ whiah loada to a 
set of nonlinaar a^uationc. Xn practiaa thaaa aquat^ona hava to ba 
solved by hill*eliMbing taehniquaa* aithar analog or dlgital» Mhiah ia 
also an ittrativa proeadura* 


Idantlf ication in a cloaad loot 


Tor a aystaii in a cloaad loot tha ramnant nit) ia no longar unaorralat« 
ad with tha ayatao input x(t), baeauaa it eirclaa around by way of 
tha faadbaek. Thia maana that tht appliaation of apan loop nothoda 
In a closad loop will load to a blaa in tha oatlmation of tho paraoatar 
Thara era two poaaibill t las to aelva thia problaot 


3.1 |rscff9mSi9D.o{.!l>t.SlSfS4.lS9B.nnt».igSS.9B.MIlir«)fOS 
se9S.i99e.!9!St!.lt9f.n«i.}} 




FTOUKt 3t 

Sy9t0m in a aloaad loop (a) and tho oquivalont 
opon loop oyotom (h). 


S07< 


If H^(p> ia tbo unblown tranafor funation» and a fcaava 

tranafor funationi than tbo oloaod loop traaafop faaotiaa N^(p> iai 



thia tranafar funation aan ba aatlaatod by ftoana of aa opoa Idop 
nothod aa illuatratod in fig* l« Tbo tranafor fimotloii b^fp) tbon 
followa from 

«i<P) ■ 

baaad on tq, (9) a nodal of N^(p) aan ba built aa a oloaod loop 
ayatan with • known tranafor funation H,(p) and an unbnoun tranofi 
funation N^(p)* In auob a ayatan a paranatar aatination taabniguo 
aan ba aaoonpliobad aooordiug to figt a or noaording to rigt I* 




ffgpbi 4t 

Foronotop oatdiNitdoii in a otoood loop 
oboro Bg(p) ia 0 bnoim apoeoaii 

Tho lattar natbodi for inatanaoi la appliad by fobannaon [a], 

Tbo natbod naotionod boforo inplloa that bnonlodgo of tbo tranafor 
funation H^(p) of tbo naobiito ia pogdirod* ^o^oooPt in nany praotio* 
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FiaUBS S: 

Mod^^ft0d vttvion cf th& m0thod cf Fi$, 4 , 

aa tittiationt this knonladga ia not avail«blo» at ioatt not to tbo 
dagrao of accuracy roquirad to gat a roUablo aatiiiato of th# 
trar jfar function M^(p)* 

In a closed loop an unblasod attiAsto can bo ebtalnod only if tbo 
remnant la sore or if tho r sonant can bs soparatod froa tho iignalo 
used for the Identification ( this can bs aebiovod by applying 
f 11 cars. rig. 6 ahows a block diagram of a aothod to obtain un« 
biaaad aatloataa. Tha input x(t) and tbs output y(t) of tbo onbnoon 
systam Hj^(p) ara flitarad in suoh a way that only thono oottponontt 
of tha algnaia x(t) and y(t) which originate in tho OMtopnnl forging 
function r(t) contrlbuts in tho parsootsr oatioation* If tbo forging 
function r<t) is a stochastic signals tha filter operation oonoioti 
of tha computation of tha eress-csaarlanca funetiona of tho foroing 
f inecion r(t) with tha input x(t> and with tha output y(t) roopootir* 
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ff 

Unbia$ 0 d popMotor §9timation for a gpoeon in a oioootf 
loop ifbori both tk§ trona/or /iinotiofio 8m(p) ond 8^(p> 
art nnktioimi 

oly. h ooro attrootivo noy of filtoring oon bo aobiovod if tbo 
ontornol foroing funotion r<t) oonoitta of a tiMi of sCnogoidgi fob 
tho gpplioition of tbo filtorg it oquivoiont to tbo ootiamtien of 
a aat of fouriar oeoffioionta. Tboroforoi «ony iovootigatora 
Cm it 7 ] bpplp • too of aingioida •• 0 foroing foootion in gnpori* 
■onto to idontify tho tronafor funation of tho hmioii oootrollor. 

In thie oloio of inroatigotiona it in important to ohoooo tbo oombor 
of ainuioido in tbo foroing funotion ouffioiontiy bigb« 00 tbot tbo 
input gppooro no 0 rondo* oigool to tho huMo oootrollor* 

* • * *.sicii<.l>sB.si!ts<.«iigi.i.<ts<niiif sli.im.ilnti 

itorting fron a difforont biobgrouadt via. tho oppliootioii of 
hinory oolti*froquoooy toot oignolo for gyoton idontifiootion (tit 
von don Boa orrivod indopondontiy at tho nothod loot indloofod. 
NoroovoPi ho olfo ohowod that it oould ho oppliod to tho gonoroliood 
modal in on opon loop opotom on wall 00 in a alotod loop eno uithsut 
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getting blat«d attlmat«8 . Hanca tha nathed juffmantioaad oan 
b« traatad a mathed for a ayatan in an apan Xorp tfitb taro ran* 
nant. This outUna will ba fallowad in tha nant darivation. 


Consldar a ayatan having a tranafar funatlon H(p)i 


H(p) s. 




In whi^h p la tha Haavlalda aparator pad/dt# Tha input xlt) ta thii 
aystam la daaerlbad by 


x( t ) s i (aj^ eeauj^tvbj^ alnuj^t) ♦ nj<t) ® x*^(t) a a^<t)( (5j 

tha output y(t) la daaerlbad byi 
n 

y(t) 5 I (Cj^ coaw|^tad|^ ainw|^t) ♦ n2<t) • y*(t> ♦ (g) 

Tha astlnatas and of tha Pouriar eelffioianta and aan 
ba obtained Promt 

9 T 2 

“ T ^ x(t) coawj^t dt 8 aj^ ♦ Y / eotOj^t dti ( 7 ) 

T T 

a y / x(t) ainwj^t ^ f / alnw^t dt (g) 

In thaaa fornulaa tha obaarvatlen time T la tha nariod of tha 
f undaTn..ital fraquaney of tha aignaia xlt) and y(t)i Similar ax* 
praasiona are valid for tha aatlmataa \ and 3^ of tha aalffialanta 
and d|^. Tha filtered algnala S*(t) and 9^(t) oan ba dafinad aat 


s I eoawj^tvbj^ alnO|^t)i 


^•(t) = 
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Hanaa tha larfar part of tha diaturbaaaaa a^(t) and ng(t) la 
filtarad out from tha sipnala x<t) and y(t)i only tha amall part 
around tha airaular fratuanaiaa ramaina* Tha ralatian Patraan tha 
dttarminiatii parta x*(t) and y*(t) af tha input n<t) and autput 
y<t) raapaevivaly aan ba daaarihad byi 

UtS^pvt.itS^p”) y^(t> • (apta^pt, »*ao^p^> n*<t*f^)i <U) 

In a madal with parnmatara d^« d^i ^l**'^n* ** 

input d*(t) a aiwilar ralatian ia titan by 

<ial^pa.#,aj^p*)9*<t) t (S^tt^pt. • •♦a^p^)t*(t*f a i(t)t <12) 


a nara aanaral way af writing this aquatian iai 

C(t) • y(t> - ■ y(t) • l^(t»f)£* (II) 

Xn Cq» (II) tha funatlon y(t) aarraapanda ta P^it) ia Bq« (12), 

Tha vaatar £ aanaista af tha aatimatas 9^ and Ij of tha uakaawn 
paranatara^ and tha taator £<t|t) aanaiata of tna nagntivta af all 
aanaltlvity funatiana 


and whiah in thin oaaa art nqual to *p^(h (t*T» and 

p^(^*(t)> raipaetitaiy, Tha eptiaal talua of t aarraapanda to tha 
aatlmatad tiaa dalay HaW| eaaaidar tha oritarian funatiaa 

if/ dt. (IH) 

0 


niaimiaatioa af thia arittrlaa funatiaa with raapaat ta tha param* 
atari to ba aatiaattd yialda a aat of aquatiaaa from whiah thaaa 
paramatara eaa ba aoltad» Pigura 7 which ia ia faat a aambiaatiaa 
af tha Piga 2 and a ahowa tha ganaral idaa of thia aathadi Tha 
mathad ia alabaratad in mera datail a# fallowa, Pafina tha aoalar 


n 



(t) dti 


<ii) 


tha vactar 
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FIGURE ?: 

VnhiaB^d pai*asi0t«i* 0$timation in a olo»$d lopp 
a gpnnratiMPd modntt 

. T 

u(T ) * f / 
an4 tha matrix 

f 

Z<t) s i / *(t,T) *^(t,T) dt. 

Than 
Z(t) * 
and 

E(a,T) s n - - u^(t)o ♦ o^ 2 (t) 0 » 

whara 

T T 

«u(t) »u<t )o , 


(U> 

( 17 ) 

( 18 ) 
( 19 ) 
( 80 ) 
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Zq, (29) can ta aoivad by an itaratieo prooadiira auoh aa tha 
Sawton'Raphson algoritho. In ordar to apply thia algeritlia» tha da* 
rl\jtive of f 2 (f> haa to ba haown« Tha darlvativa of f 2 <t) ia givaa 


— 7 ; — » |z(t)| (2|^»u^(t)) Adi 

2 

* 8(t)i^ u<t)» 

2 2 

♦ U^(t)2-- <Adl Z(T))JI<T)] - ^(|Z(T)|}j(^<t> Ad} 2<t)£<t). 


Cq. (29) will hava an infinlta nuobar of aolutionat but 
ow ing to the fact that In praetioa tha range of poaalbla aolutiona 
for the value of t can ba given t nornally an unanbiguoua aolutioo 
can be found* The starting point for tha iteration proaadura oan ba 
chosen by eolving £q, (24) for a anall nuabar of value# of t in 
the range of interest » so that Eg. (29) oan ba solved, finally » 
the paraneter vector a can be coaputad fron Eg. (23). 

For a given aodel and for known aatinatas of tha fouriar coaffioiants 
of the signals K(t) and y(t)| all elanents of tha vaotora and aatriaas 
mentioned in Ege (29) and Oo) nan ba conputad. Per aKanpla» ift 

then Eg. (12) can be written ass 

t(t) » f*(t) - {8^s"(t-?,) t 3jP(a*(t-?,» - („) 

wcich means that ia this ease the alsnants C^(t«t) of tha vaster 


c.(t,T) a S^(t-T) a [ (I. eosa.(fT) ♦ S,^ ain«.(t*r>)| 


n ^ 

Cj(tpt) s p{R^(t-t)> a J ees«j^(t*T) - sin«j^(T«T )) ( (34) 
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Ih 

s -p(9*(t)> S 

Tha alamants of tha vaster ^(t) and tha natrix E(t) aan ba ealau- 
latad in a way ainilar to tha ana given balawt 
T 

t f ^ dt • 

t - 4 p(9*«t)), 

0 

»!»<»> • * ••••k* ♦ <Mk*Mk> •»"«k»>- 

frofli tha foregoing it follow# that all infornatian naaasaary for 
tha aatimatian of tha par«»#t#r# S# available if tha aatinatas of 
tha faupiar aoaffiaiant# of tha signals involved are known. 


for systans having more than one in;:* and nara than ana output t 
tha idantifieatien of unknown transfer funations is lass straighf 
forward. Con#idar» for instanast the systan of Pig. 8 whiah ra- 
prasant# a nan*bieyala systan whars tha rider baa to parfarn two 
tasks vis. tha stabilisation of the bioyalo and the following of a 
given traek. Xn thia tyatan inpute and outputs of tha unknawr trans- 
fer funationa H^(v)i HjCv)! H^Cv) and H|^(v) are oaupladt hat only 
within tha hunan operator* but aiao within the bioyclo ainulatar. 

Tha aignala R^(v) and are introduoad axtaraally aa feroing 

funetiena. If tha ralatlana batwaan inputa and outputs of tha bicycle 
iinulatar are linear* tban all signals in tbs systsn asn bs dss* 
oribsd ss linssr funstions of tbs feroing functions R^(v'.* R^(v) and the 
rsnnsnts b^(v) snd N^Cv). far instanaat 

W<v) i P^(v)Rj(w> ♦ P2<v)k2<v) ♦ Pj(v)h^{v) ♦ P^(v>h2<v>» ^37) 

X(v) e Q^(v)R^(v) V 02 (v)R 2 <v) t v Q^(v)H 2 (v)* (31) 

where the tranafar functlo Pj^Cv)* p 2 <v)* p 3 <v)* F^(v>* Q^(v)* O^Cv), 


-516 
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PIGVRB 8: 

Bxampla o/ a Ttultitocp 

G (v) and G^(v) describe the relations between the external input# 

R (v) and on the one hand and the signals W<v) and X(v) on the 

other, rro* Fig. 8 and Cqs ( 37 ) and (ag) it follow# thati 

y(v) r (Kj(v)r^(v)tH2(v)6^(v))R^(v) ♦ {H^< v)rj(u )tHj( u)Oj<u))Rj(v) 

T UtH^(v)F3(v)tHj(v)G3(M)}H^<v) ♦ {H^( v)f^( v )eH2< v) 0 j^( v)}»3< v ) | 

(39) 

2(v) s (Hjl j)F^<v)tH^(v)6^(v))R3(w) ♦ (H3(u)Fj< v)eH^(g)«2<w))Rj< v> 

t (H 3 (v)F 3 (v>tH^(v) 03 <v))ll 3 (w) ♦ (ieH3rv)P^(u;tB^<v)0j^(w))Hj( v)t 

(40) 

Suppose that for each of the signals W(v)i X(v)» Tf(v) and 2 (w) the 
conponents W^(v>* X^(v)» X^Cv)! Y^(v)* 2 ^(v) and Z^Cv)! 

originating In the external test signals R^iv) and Rjtv) can bS 
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id* Than a set of filtered eignale oan be dietinguiebed« vis*i 


Wjlv) 8 F3 (v)R3(v)| 

Wj(v) 8 Fj(v)R3{v)j 
X^(v) 8 03 <v)R^<v)* 

XjCv) a 83 (v)Rj(v)j 

y^Cv) 8 (H3 <v)Fj^(v) e H 3 (v) 03 (v)>R^(v)i 
VjCv) 8 <H3(v)f3(v) ♦ H 3 (v) 03 (v))R 2 (v)t 
Z^(v) 8 (HjCvlF^fv) t H^(v)C3<v)>R^(v)t 
Zj(v) 8 (HjCvlFjCv) e H^(v) 03 (v))R 3 (u). 


By eliminating F^Cv), *** 

can be redueed to a set of a eguatiene from whieb the 4 naknotm 
transfer funetiene H 3 <v) end B^(v) oan be aelved* Ae 

an example I 

X*(v)T,(v)»X,(v)y,(v) 




i}<4-«;Nr78T 


Hj(w) » 


»,(v) 

y '(v) 

^ OTTT* 


Separation of the eemponents of the signals originating in the two 
forcing functions is possible when both these tost signals ais 
oempesed of a nanber of tinttseida* Xn order to distinguish between 
the eemponents originating in eaoh of the two test aignale it is 
neosseary that no eemmen fregueneles esour in both teat aignala* 
However^ new the problem arises that applieatlen of Eqa (42) or 
(43) is not possible booause for a given frequenoy either the 
signals with the index 1 or those with the index 3 or both are 
sero. 2n general the transfer funetiens eenoldored here are euffi- 
oiently smooth » itO. the tranafer fnnstiena oan be ooneidorod to 
be oenetant within a frequenoy ranga Av * 1/T* Row the problem oan 
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b* sol»*d «» indicated balow. ,iBU*olda with 

Cboo.a aft th. .i.«.l H,(v) which 



Than Lt it la asauaad thati 

I, »• 

i « tran«f«r function* Ha(vi* **2' ”3' 

and H^(\), It i „ ,1 \ 1 - Ea (as)* Thaaa daocwplad IB- 

output. Xlk. "i’*''’ ia\h. para-tar a.tl.atloB -thcd 

puta and outputs c -.ibility to obtain tha uaeoupxad 

,..u.a t. .. wpiT ~ 

.» .. — >• - 

et 8 l [10] . 

u .»..w rr:::” ” 

b* b*el..t.d P. ,, bbbbtPb.. A paA.pAp.pp 

;:;;:n.:-r‘T..';.:r..T.r:rp r—. ...... 

.... .( A01..1... ... P. “• 
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